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ARTICLE INFO ABSTRACT 



Multiaxial fatigue is a very important physical phenomenon in several mechanical components. Fatigue 
life study under cyclic stresses is of utmost importance to avoid unexpected failure of equipments, vehi- 
cles or structures. Among several fatigue characterization tools, a correct definition of a loading cycle 
under multiaxial fatigue loading conditions shows to be crucial to estimate multiaxial fatigue life. 

The aim of this work is to achieve a correct definition for a multiaxial fatigue loading cycle and accom- 
plish a multiaxial fatigue model to estimate block's fatigue life under multiaxial loading conditions. To 
reach this goal, several loading paths were carried out using the 42CrMo4 low alloy steel under different 
loading conditions. Sequential, proportional, non-proportional and asynchronous loading effects were 
modulated through eleven loading blocks. Furthermore, two models were proposed: a cycle counting 
method and a fatigue life evaluation criterion. The results from the proposed models were correlated with 
the fatigue data and compared with two well known cycle counting models: the Bannantine and Socie 
and the Wang and Brown criteria. The proposed models were successfully validated by experimental 
data. Results show that the new proposals lead to an improved multiaxial fatigue characterization under 
complex loading conditions. 

© 2014 Elsevier Ltd. All rights reserved. 
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1. Introduction 

In simple fatigue loading cases the cycle definition using an 
equivalent stress is very prompt and the inherent fatigue life cor- 
relations provides acceptable results. However, in complex loading 
histories, the fatigue life estimations under equivalent stress ap- 
proaches, most of the times, are not in agreement with experimen- 
tal data. Equivalent stress approaches are unsuitable to 
characterize directly the block's fatigue damage, because they yield 
non-conservative fatigue live estimations. In other words, the load- 
ing block's fatigue damage is greater than the one estimated by the 
greatest equivalent stress found in that same loading block. 

Therefore, block's fatigue damage characterization must con- 
sider what happens during the entire block's loading history. Mul- 
tiaxial cumulative damage is commonly calculated using the 
Miner's rule, or its different versions, as is used in uniaxial loading 
conditions. 

In addition, Miner's rule considers a damage parameter in asso- 
ciation with a cycle counting method [1 ]. Damage parameters must 
capture the physical fatigue behavior to allow considering the 
cumulative damage as a cycle damage summation [2-5]. 
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Equivalent damage parameters, such as equivalent stress, keep 
the time relation between multiaxial load components, which is 
a plus, but they have some shortcomings. For instance, time histo- 
ries of an equivalent stress are always positive and due to that, 
their time history is unsuitable to capture the whole physical dam- 
age process. 

The load sign of multiaxial load components is essential infor- 
mation because it identifies a different damage behavior. For in- 
stance, compression and tension stress states have negative and 
positive sign and distinct fatigue damage behaviors. Compression 
stresses tend to slow down the fatigue damage and the tension 
ones tend to speed up the damage process. Moreover, SN damage 
curves are established based on a pair of positive and negative 
reversals i.e. a tension and compression load pair. Each uniaxial 
reversal contributes in a different way to the SN unitary damage 
composed by the compression/tension damage pair, in the case 
of an axial SN curve. 

Therefore, a rainflow cycle counting over the time variation of 
an equivalent stress does not capture whole the fatigue damage 
behavior [6,7], 

Nowadays, new approaches on cycle counting are very few [5- 
10]. Well-known cycle counting methods use the rainflow method 
or its variants. They have acceptable results under uniaxial load- 
ings, but under multiaxial loading conditions the fatigue life esti- 
mations yields poor results. Some cycle counting criteria are 
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critical plane based methods, they use the rainflow method applied 
to the stress/strain time variation, on each projection plane, to ac- 
count loading reversals. For instance, Bannantine and Socie (BS) 
proposed a cycle counting criterion that accounts the number of 
loading reversals, in the critical plane, through a rainflow routine 
[11]. 

Subsequently, a damage parameter is determined, at each 
reversal, using a critical plane model according to the axial or shear 
based damage approaches. BS method has some shortcomings in 
damage characterization, because the contribution of the axial 
and shear loading components to the fatigue damage are discon- 
nectedly computed. Therefore, the time relation between shear 
and axial load components is ignored, and the joint damage effect 
is not captured. Wang and Brown (WB) [2] proposed a new cycle 
counting method to overcome the sign lost shortcoming in their 
equivalent strain approach. WB cycle counting criterion is based 
in the von Mises relative strains and in the classic rainflow cycle 
counting method. 

The loading cycles are extracted from the von Mises relative 
strain time variation and the WB damage parameter is computed 
for every extracted cycle. Then, an accumulated criterion must be 
used to calculate the loading block's total damage in order to esti- 
mate block's fatigue life. The WB criterion has one main shortcom- 
ing, which is the possibility to miss the greatest stress range during 
the loading time history [12,13]. 

Meggiolaro and de Castro, proposed one modification to the 
Wang and Brown criterion to overcome this shortcoming. The 
change was to expand the WB criterion to the five-dimensional 
Euclidean space to ensure that the highest relative strain ampli- 
tude is always accounted [6,7]. Wei and Dong proposed an equiv- 
alent stress cycle counting method using the time evolution of the 
equivalent stress in the von Mises stress space. The equivalent 
stress is mapped to extract loading cycles and overcome the miss- 
ing loading path-dependency verified on the WB criterion [12-14]. 

In this paper is studied eleven loading blocks with several dif- 
ferent loading effects using the Bannantine and Socie, Wang and 
Brown and a new proposed multiaxial cycle counting method. 
The new cycle counting method, the Stress Scale Factor (SSF) vir- 
tual cycle counting, is based on the SSF equivalent shear stress 
early proposed by the authors in [15]. 

Moreover, a new criterion to evaluate fatigue life under com- 
plex loadings is proposed, this method uses the new cycle counting 
method proposed here and is an update of the SSF equivalent shear 
stress to block loading conditions. 

In order to validate the proposed methods it was performed a 
fatigue life correlation using fatigue data from eleven loading 
blocks. Moreover, the proposed models estimations were com- 
pared with two state-of-the-art models estimations in this field. 
Results show that the proposed cycle counting and fatigue life esti- 
mation approaches yields better results. In addition, the proposed 
models computational work is far less than the one necessary in BS 
and WB fatigue life evaluation. 

2. Theoretical development 

To estimate fatigue lives from loading blocks is necessary to 
consider a multidisciplinary approach in fatigue damage character- 
ization. Multiaxial block loadings have much more intricate load- 
ing histories than the reference ones used to set up SN curves. To 
achieve block loading damage is necessary to enter with three 
main fatigue approaches: a cycle counting method, a damage crite- 
rion and a damage accumulation model. The damage criterion is 
the base stone to set up random fatigue. The damage parameter 
must capture the fatigue damage behavior to allow set up a cycle 
counting method and an accumulation model. 



2. J. Multiaxial cycle counting 

Uniaxial cycle counting is a well-understood matter. Hysteresis 
analysis to find a load cycle is a method that has proved to be reli- 
able in uniaxial fatigue damage characterization. However, under 
multiaxial loading conditions such prompt method does not exist. 
Efforts has been done to adapt the uniaxial cycle counting method 
to the multiaxial loading histories [16-19]. BS and WB cycle count- 
ing methods are an example of that; they use the rainflow para- 
digm to extract the loading cycles within a multiaxial loading 
history. 

2.2. Bannantine and Socie (BS) 

Bannantine and Socie [2] proposed a multiaxial cumulative 
damage criterion based on a mix of critical plane damage parame- 
ter, rainflow cycle counting method and Miner's rule. The main 
concept behind BS criterion is based on the experimental evidence 
that some materials are more sensitive to axial strains than to 
shear strains and vice versa. 

Therefore, authors concluded that the axial and shear strains 
have different damage scales in the same material. Bannantine 
and Socie considered that the block damage under multiaxial load- 
ing conditions could be estimated through axial strains or shear 
strains. 

Thus, the criterion of choice is based on the material sensitivity 
to axial or shear damage. To evaluate complex loadings Bannantine 
and Socie proposed that the multiaxial axial and shear strains com- 
ponents could be treated as uniaxial loads in the critical plane 
search. Therefore, the multiaxial loading is projected into the can- 
didate plane, and then the in-plane and normal strains are sepa- 
rately computed with a rainflow method. 

After that the axial and shear accumulated damage is calculated 
for each plane. The plane with the highest accumulated damage is 
the critical plane. BS approach states that the rainflow cycle count- 
ing on each plane can be associated with the hysteresis loops in 
that same plane. The BS approach uses the original main paradigm 
behind the rainflow cycle counting methodology [11,20-22] to 
capture the critical plane damage. The damage accumulation at 
each plane is accounted through any cumulative damage model. 
The most used is the Miner's rule, where is performed by a linear 
summation of each loading cycle damage. Several accumulative 
damage approaches based on the Miner's rule can be found in lit- 
erature [1,23]; they are mainly a non-linear versions of the Miner's 
original law. In the authors' opinion, the non-linearity in damage 
accumulation must be left to the damage parameter since SN 
curves are already non-linear. Therefore, in this study the linear 
Miner's rule is adopted. The BS damage parameter is determined 
based on a specific critical plane criterion. For instance, if the axial 
strains are cycle counted, then a critical plane model based on nor- 
mal strains must be considered. Likewise, if the in-plane strains 
(shear strains) are cycle counted then a critical plane based in 
shear strains must be considered. For example, reversals damage 
can be accounted using Eq. (1 ) in shear strains and in axial strains 
can be used Eq. (2). 

Fatemi-Socie (F- Socie) damage parameter: 

"f{% • (l + k ■ ^) } = |(2N ; )^ + y' f (2N f r (1) 

where A)' max /2 is the maximum shear strain amplitude on a 0 plane, 
ffn.max is the maximum normal stress on that plane, a y is the mate- 
rial monotonic yield strength and k is a material constant, k = 1 .0 in 
this case. G is the Shear modulus, 2/Y/is the fatigue life and t f , b, y' f , y 
and c are material's constants (cyclic properties). 
SWT damage parameter: 
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= -^(2N f ) 2b + rj'j 2 s'(2N f ) b+c 



(2) 



where a n is the normal stress on a plane 0, and Aei is the principal 
strain range on that plane. £ is the Young's modulus, 2N f is the fati- 
gue life and %'*, b, e' and c are material's constants (cyclic properties). 
Socie and Marquis [2], suggests that the damage within a block 
loading can be estimated by adding in each candidate plane the 
damage from normal and shear strains. The main idea is to get 
the total fatigue damage through a combination of two different 
damage behaviors. Eq. (3) presents this idea, where the largest fati- 
gue damage occurs in the plane with the smallest value of N(0), (i.e. 
the critical plane). 



1 



1 



1 



N(fl) N a (ff) N t (0) 



(3) 



N a (0) and N T (8) are the fatigue life estimations obtained through a 
critical plane based on axial and shear strains, respectively. The Socie 
and Marquis approach aforementioned has one main shortcoming 
because it considers twice the same axial damage contribution. In 
fact, if the total fatigue damage is accounted based on Fatemi-Socie 
and SWT damage criteria, one fraction of axial damage is double 
counted. This is so because F-Socie criterion corrects the shear strain 
amplitude with the greatest axial stress. Usually, the critical plane cri- 
teria adopt the highest damage parameter to estimate fatigue lives. 
However, there are differences between the highest damage and 
the highest accumulated damage in fatigue life estimations. The high- 
est damage parameter considers just the load reversal with the high- 
est amplitude within a loading block, at a specific plane. Therefore, 
damage from other reversals is not accounted, but in accumulated 
damage criteria every reversal in a specific plane is accounted. In or- 
der to check the BS criterion in this study, the block damage is ac- 
counted by identifying the critical plane with the highest 
accumulated damage based on shear or normal strains. Moreover, 
was adopted here the ASTM rainflow cycle counting criterion and 
the Miner's rule to achieve the total fatigue damage at a specific plane. 

Fig. 1 shows the Bannatine and Socie method to estimate a 
block fatigue life. 



2.3. Wang and Brown 

Wang and Brown [24] proposed a multiaxial cycle counting 
model where reversals within a loading block are extracted 
through the von Mises equivalent strain/stress. This approach just 
identifies the region from the loading path that belongs to each 
reversal. 

A damage parameter can be applied to each one of the reversals 
identified as if it were a new and unique loading path. The von 
Mises equivalent strain/stress has attached the loading sign loss 
shortcoming i.e. the equivalent strain/stress is always positive 
regardless compression, tension or shear directions. 

This is an important matter because completely reversed load is 
transformed into a zero to maximum creating a virtual mean stress 
on the loading path that does not exist. 

In order to avoid this shortcoming Wang and Brown have pro- 
posed to make a rainflow cycle counting to the relative von Mises 
strains. To extract the reversals within a loading block, the multi- 
axial loading history is calculated to get the von Mises equivalent 
strain time evolution. After that, the maximum equivalent strain 
location is identified to move the time origin to that same location. 
Identifying the maximum equivalent strain location is the preli- 
minary modus operandi to carry out the Wang and Brown cycle 
counting method. The time evolution of the relative equivalent 
strain is then calculated by computing the norm between the 
instantaneous equivalent strain and the maximum equivalent 
strain value. After that, a rainflow cycle counting is performed into 
the relative equivalent strain time evolution. Later, the reversal is 
extracted and the relative equivalent strain history is updated by 
removing the information inherent to that same cycle. The remain- 
ing loading history is computed in the same way until finish all 
available loading branches. 

A map is obtained with the WB method; it splits the loading 
path into loading regions i.e. loading reversals. For each reversal 
a critical plane search is performed to find out the damage param- 
eter inherent to that cycle. 

Eq. (4) shows the Brown Miller critical plane criterion typically 
used with the Wang and Brown cycle counting method. 



Axial and shear stress 
plane projections 




Axial strain 
rainflow cycle 
extraction 




Shear strain 
rainflow cycle 
extraction 




Rainflow axial strain results 



SWT parameter at each reversal 



Nf for each SWT parameter 



Miner's rule at each plane 
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Rainflow shear strain results 



F-S parameter at each reversal 
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/.each reversal 



Block fatigue damage 
is given by: 



D e =Max( D e , D 0 ) 

axial shear 



Block fatigue life 
estimation is given by: 



N 



f .Block 



Fig. 1. Bannantine and Socie methodology. 
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l+v' + (l+v')S 



2<7„, r 



-(2N / ) b + e}(2N / ) c 



(4) 



where Ay max and 5e„ are the shear strain range and the normal 
strain excursion between two turning points of shear strain com- 
puted on the maximum shear strain plane. The S parameter near fa- 
tigue limit is equal to 0.3, for half fatigue life S takes the value 0.7, 
[2]. 

£ is the Young's modulus, o n , mean is the mean normal stress 
within the loading path, 2Nf is the fatigue life and x',, b, e' and c 
are material's constants (cyclic properties). 

Subsequently, a cumulative damage criterion must be used to 
account the block damage to estimate the inherent fatigue life 
[24]. Fig. 2 shows an example of the Wang and Brown method 



where it is found out the reversals within a loading block. Fig. 2a 
shows the axial and shear loading components time evolution. 
Afterwards the von Mises strain/strain time evolution is computed 
see Fig. 2b. Next step is to find the first von Mises strain/strain 
maximum value and translate the graph origin to that point. In 
Fig. 2c is shown the referential translation for the first maximum 
von Mises equivalent stress named as point A. 

At this point the relative von Mises time history ((T iretative ) is 
determined by subtracting, in the translated loading path, the axial 
and shear components at origin. Eq. (5) shows the numeric method 
to get the relative time history. As a result, the relative von Mises 
time histories always start at zero. 



<7i,retotf«*(t) = V - ff rmx) 2 + 3 ■ (T,(t) - T nux ) 2 



(5) 
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where cr, and x\ are the instantaneous axial and shear stresses, 
respectively. o" max and T max are the maximum axial and shear stres- 
ses found within the loading histoiy. 

To extract a reversal it is selected on the relative von Mises time 
histoiy the loading path branch that starts to increase from zero. 
When the relative time history starts to decrease, the reversal 
extraction goes on the horizontal until the time history starts again 
to increase. If the horizontal extraction does not find any increasing 
relative time history branch, then the reversal ends at the point 
where the extraction initiates to become horizontal. 

Fig. 2d shows the 1- reversal extraction represented by line AB; 
at this stage, no highest point than point B is found, so the reversal 
was extracted. Later, the axial and shear components inherent to 
the reversal extracted must be removed from the loading block 
path. Then the procedure starts again with the updated axial and 
shear stress loading path; just one reversal is extracted for each 
one relative time history. In Fig. 2e is shown that point B is the 
new relative minimum. Continuing the reversal extraction process 
it is achieved point C; at this point the von Mises time history de- 
creases, then the dot line goes straight to find the point C 

The reversal extraction continues to the time histoiy end. The 
second reversal was extracted with the points BCC'D. Then the ax- 
ial and shear stress loading paths are updated again. The new von 
Mises relative minimum is located at point C from where grows 
until it reaches point E, being extracted the 3rd reversal, CE, see 
Fig. 2f. 

Next, the axial and shear stress is updated by removing the 
reversal CE axial and shear components from the original loading 
path. After that the new relative von Mises time evolution is calcu- 
lated, see Fig. 2g. 

The reversal extraction starts at point E and increases to reach 
point F. At this stage, entire loading path was updated and it was 
found eveiy reversal within the loading block. In this example, 
the loading block time history was separated in four sections, i.e. 
four reversals. 

2.3.1. Wang and Brown block fatigue life estimation 

Each one of the achieved sub-loading paths (reversal zones) 
must be considered as independent loading paths. Then, for each 
reversal, the damage criterion must be calculated and the inherent 
fatigue life estimated. 

Afterwards, Miner's rule is used to compute the block accumu- 
lated damage and the block fatigue life can be estimated, using 
Eq. (6): 



D h 



^reversals 

- E 



1 



•fJ 



(6) 



where Nfj is the reversal fatigue life estimation and D b i ock is the to- 
tal damage within a loading block. Thus, the block fatigue life can be 
estimated as follows, Eq. (7): 



1 



(7) 



instantaneous damage contributions of axial and shear loadings, 
which is a huge advantage. 



2.4.1. SSF equivalent stress approach 

The SSF equivalent stress approach was carried out by the 
authors [15 j. The proposed approach considers that the stress 
amplitude ratio and the stress loading level have a huge influence 
in the material fatigue strength. These influences were accounted 
through the SSF function that transforms an axial damage into a 
shear one. With the SSF equivalent shear stress it is possible to esti- 
mate fatigue lives using Eq. (8). 



max(T + ssf ■ a) = A(NA 

block 



(8) 



where % and a are the instantaneous load components, A and b are 
the power law trend line components of the torsional fatigue life re- 
sults (torsion SN curve), and N/ is the fatigue strength under certain 
stress level. The SSF function is given in Eq. (9): 



ssf (a, X) = a + b ■ a + c ■ a 2 + d ■ a 3 +f ■ ) 2 + g ■ '/? + h- i 4 + i 

• ;. 5 



(9) 



where a is the axial component of the biaxial loading and X = x/a is 
the stress amplitude ratio. Constants from "a" to "i" are determined 
through experimental tests. The SSF function is a material fatigue 
property and must be determined by experimental tests [15]. 



2.4.2. SSF virtual cycle counting 

Cycle counting methods are used in fatigue life evaluation be- 
cause the maximum damage parameter within a loading block 
does not capture the block's fatigue damage. Therefore, block's fa- 
tigue life estimations obtained solely with a maximum damage 
parameter yields poor results. Here it is presented the SSF virtual 
cycle counting method. In this new approach are not analyzed hys- 
teresis loops and due to that it was named as virtual cycle count- 
ing. However, the new method is physically based. It relates the 
maximum damage parameter, within a loading block, with the 
block's total damage. 

Fig. 3 shows the SSF virtual cycle counting method in associa- 
tion with the block's fatigue life estimation method. Fig. 3a pre- 
sents the axial and shear stress components time variation of a 
multiaxial loading which are computed to get the SSF time histoiy, 
see Fig. 3b. Next, in Fig. 3c is selected, from the SSF time histoiy, 
the highest SSF value and the SSF stress value at every peak and 
valley found between two consecutive SSF zero stress points. 

The highest SSF value is the block's damage reference to achieve 
the block's damage, see Fig. 3d and e. The virtual cycle counting is 
achieved by adding each SSF peak/valley absolute values and divid- 
ing the result by the block's damage reference. Eq. (10) shows the 
virtual cycle counting expression. 



where N Wock is the block fatigue life estimation. 



vcc -. 



Eabs(x] 



peak, valley 



2 ■ T m axMock 



(10) 



2.4. New cycle counting and block's fatigue life evaluation criterion 

Authors propose here a new criterion to estimate loading 
block's fatigue life. Moreover, a new cycle counting method to ac- 
count loading block's fatigue damage is also proposed, i.e. a non- 
rainflow cycle counting method based on the SSF equivalent shear 
stress is presented. In the author's opinion, equivalent stress ap- 
proaches are suitable to account block's fatigue damage. Equiva- 
lent stress criteria packed in one damage parameter the 



where vcc is the virtual cycle count inherent to a loading block and x 
is the SSF equivalent shear stress at each peak/valley. When the 
loading block is defined by two SSF reversals, i.e. one peak and 
one valley, the vcc is equal to 1. In this way SSF virtual cycle count- 
ing accounts for 1 cycle loading cases, that usually are the reference 
ones. 

Carrying out the proposed multiaxial cycle counting method it 
is much easier than other cycle counting methods such as BS and 
WB. 
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2.4.3. New block's fatigue life criterion 

The block's fatigue life estimation is achieved by using the 
block's reference damage and the virtual cycle count as follows, 
Eq.(ll): 

N/^ = %= (ID 

where Nf_ b i ock is the block's fatigue life estimation and N x max is the 
fatigue life estimation reference damage determined by Eq. (12): 

% t max = (^±y (12) 

where A and b, are the power law regression components obtained 
from pure torsion fatigue data. 



3. Materials and methods 

3.1. Material 

A 42CrMo4 quenched and tempered high strength steel was 
used in this work. The test specimens were machined from rods 
with a 25 mm diameter. 

The specimen's geometric shape was a solid hourglass with a 
6.35 mm diameter at middle length [15]. 

Specimens were inspected and manually polished using sand- 
papers of decreasing grain size, from P200 to PI 200. Fatigue tests 
were carried out through a multiaxial servo-hydraulic machine un- 
der load control at room temperature. 

Tests were performed under axial and torsion fully reversed 
conditions as shown in Fig. 4. Testing frequency was 5 Hz. 



3.2. Loading blocks 

To do the comparative study between the accumulated damage 
criteria presented in the earlier sections, eleven loading blocks 
were tested and analyzed. 

Fig. 4 shows the loading blocks on the axial and shear stress 
time history evolution. 

4. Results and discussion 

4.3. S-N results 

Table 1 shows the experimental fatigue life results for the se- 
lected loading blocks. The Axial and Shear labels, presented in Ta- 
ble 1, are the highest stress in axial and shear loading channels, 
respectively. Nf is the number of cycles to failure and (ro) is a 
run out test. 

4.2. Theoretical cycle counting results 

Theoretical cycle counting results, for the selected criteria, are 
presented in Figs. 5-1 5. For each loading block, the selected criteria 
and the new approach results are presented in the same figure to 
turn easier the comparative analysis. 

4.2.2. How to interpret the results 

Cycle counting results are presented here in different ways be- 
cause each criterion have distinct approaches and different ways to 
yield the results. However, the number of cycles within each load- 
ing block is focused here as a common result between the selected 
criteria. In Figs. 5-15, the sub-figures (a-c) are related to each 
other. The loading path in sub-figure (a) is presented in the von 
Mises stress space, from which is analyzed the critical plane 
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(i) 



(j) 




Fig. 4. Axial and shear stress time evolution of the selected loading blocks. 



estimations and the loading path shape. In experiments, specimens 
test had their longitudinal direction lined up with the a axis, 
shown in sub-figures (a), and 0° crack direction is parallel to z axis. 
Bannantine and Socie results are shown in sub-figures (b and c). To 
estimate the loading block fatigue life the accumulated damage 
was determined and represented by a solid line at each loaded 
plane. Furthermore, a dotted line represents the BS reversals num- 
ber at each loaded plane. Moreover, sub-figures (d and e) are re- 
lated to each other. Sub-figure (d) shows the axial and shear 
stresses time variation and sub-figure (e) shows the Wang and 
Brown criterion results. WB reversals number is a time dependent 
result and the most damaging reversal occurrence can be identified 
in the stress-time loading history presented in sub-figure (d). Like- 
wise, the SSF equivalent stress time variation, shown in sub-figure 
(f) can be compared with sub-figure (d) where it is shown the load- 
ing block stress-time variation. Furthermore, the critical plane 
reversals number shown in sub-figures (b and c) have a different 
paradigm than the WB block reversal number presented in sub-fig- 
ure (e). In critical plane approach, the reversals results are related 
to the number of reversals computed in each plane. On the other 
hand, Wang and Brown reversal results show the reversal time in 
the loading block period. 

4.2.2. Results analysis 

Fig. 5 shows the results of the loading Case 1. Case 1 is a sequen- 
tial loading where a sinusoidal axial load is followed by a 



sinusoidal shear load s/3 times less than the axial one. Regarding 
BS results, the SWT and F-Socie criteria estimate the critical plane 
direction at 0° with 4 and 3 reversals, respectively. The experimen- 
tal critical plane for loading Case 1 is -5°. Considering the critical 
plane estimation calculated without accumulated damage 
approach, the SWT and F-Socie critical plane estimations were 0° 
in SWT and ±90°; 0° in F-Socie. In Fig. 5b, at -45° is shown that 
the loaded plane has 6 reversals, however the inherent damage 
parameter is near zero. 

Moreover, the maximum damage parameter was found at 0° 
where the reversals number is lower than the one found at -45°. 
It shows that the number of reversals must be always kept up in 
association with the damage parameter. Regarding WB results it 
can be concluded that Case 1 has 4 reversals. The most loaded 
reversal is the 2nd one which is divided into two time periods, 
see Fig. 5e. Comparing the WB results with the block stress-time 
history, see Fig. 5d, it can be concluded that the most damaging 
period is between 1 6 and 33 s and from 57 to 65 s. Therefore, 
WB most loaded reversal, 2nd one, outcome from the axial com- 
pressive and from the last shear loading branch. SSF time variation 
presented in Fig. 5f is very similar to the stress-time variation pre- 
sented in Fig. 5d having the zero stress points at the same time. The 
SSF cycle counting method extracts 3.4 reversals and the rainflow 
method extracts from the SSF time variation 5 reversals. 

Fig. 6 presents the loading Case 2 results. Case 2 is a sequential 
loading block like the loading Case 1 with the particularity that the 
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Table 1 

42CrMo4 fatigue life results. 

Case Axial Shear Nj 



shear stress amplitude is now greater than the axial one. The SWT 
critical plane estimation is ±45° where the BS criterion yields 6 and 
5 reversals, respectively reaching the same accumulated damage in 
each critical plane. F-Socie criterion estimates the critical plane at 
0° from which BS yields 3 reversals. 



The experimental critical plane for loading Case 2 is 0°. Consid- 
ering the critical plane estimation calculated without accumulated 
damage approach, the SWT and F-Socie critical plane estimations 
were ±45° in SWT and ±90°; 0° in F-Socie. WB results are very alike 
to the ones verified in the loading Case 1 , where the greatest dam- 
age occurs at the 2nd reversal. 

The SSF time variation and the block stress-time variation have 
both zero stress points at the same time. SSF criterion extracts 3.6 
reversals from Case 2 and rainflow criterion applied to the SSF time 
evolution extracts 5 reversals. 

Fig. 7 shows the results for the loading Case 3, which is com- 
posed of several proportional loading branches, with different 
sequential loading paths, see Fig. 7a. BS extracts 7 reversals in both 
SWT and F-Socie critical plane criteria. SWT results show that the 
greatest accumulated damage occurs at ±30°, moreover for the F- 
Socie criterion the critical plane estimation is ±1 6°. The experimen- 
tal critical plane for loading Case 3 is 29°. 

Considering the critical plane estimation and neglecting the 
accumulated damage approach, the SWT and F-Socie critical plane 
estimations are ±25° in SWT and ±25°; ±69° in F-Socie. 

WB results identify the most damaging reversal at the last load- 
ing branch and extracted 6 reversals from the loading block. 

SSF accumulated damage criterion yields 4.2 reversals. The rain- 
flow results yields 5 reversals for the SSF time variation. 

Fig. 8 presents the results of loading Case 4. This loading block is 
similar to the loading Case 3 but with different loading sequence. 
Case 4 does not have a fully reversed stress-time variation likewise 
Case 3, see sub-figure a. BS results are presented in Figs. 8b and c. 
SWT and F-Socie yields the same reversals amount and the estima- 
tions to the crack initiation direction were ±29° in SWT and ±17 in 
F-Socie. Experimental results of loading Case 4 show the crack ini- 
tiation plane occurs at 31°. The critical plane estimations without 
accumulative damage approach were ±25° in SWT and ±21°/±69° 
in F-Socie. From the results it can be concluded that the loading se- 
quence change between Cases 3 and 4 do not affect the critical 
plane estimation and experimental crack initiation direction. WB 
yields 6 reversals, see Fig. 8e, which is the same value extracted 
by the BS method. WB results suggest the last reversal as the most 
damaging one. SSF extracts 4 reversals and the rainflow applied to 
the SSF time change yields 6 reversals. 

Case 5 results are summarized in Fig. 9. This loading block is 
composed of several proportional loading branches with different 
stress amplitude ratios (SAR). The sequence number presented in 
the loading path gives the block's load sequence; see Fig. 9a. BS 
yields 10 reversals under SWT and F-Socie damage criteria. 

The critical plane direction estimated under the BS accumulated 
damage approach was 4° for SWT and 0° for F-Socie. Experimental 
results show that crack initiation plane occurs at 0°. Critical plane re- 
sults obtained without a damage accumulation approach estimates 
the critical plane direction at 0° in SWT and F-Socie criteria. WB 
method extracts 9 reversals from loading Case 5, minus 1 reversal 
than the ones obtained by the critical plane criteria in association 
with the BS approach. SSF criterion extracts 6.8 reversals and the 
rainflow method extracts 10 reversals from the SSF time variation. 

Fig. 10 summarizes the cycle counting results of loading Case 6. 
Case 6 is similar to Case 5, the difference between them is related 
to the load sequence presented in Fig. 10a. At loading Case 6, BS 
yields 10 reversals for SWT and 11 reversals for F-Socie at the most 
loaded plane (higher damage parameter). It can be concluded from 
Fig. 10b and c that the critical plane is not the one that has the 
greatest amount of reversals. 

Estimations for the critical plane directions with BS accumu- 
lated damage criterion were -2° for SWT and 0° for F-Socie. The 
experimental critical plane direction was -22°. 

Critical plane estimations without using a damage accumula- 
tion approach were -25°; 23° in SWT and -23°; 20° F-Socie 
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Fig. 6. Case 2: (a) von Mises stress space, (b) BS results with SWT, (c) BS results with F-Socie, (d) Block stress-time change, (e) WB results and (f) SSF time change. 



criteria. WB method suggests the reversal number 8 as the most 
damaging one and extracts 14 reversals from the loading block. 
SSF method extracts 6.8 reversals and the rainflow method applied 
to the SSF time evolution yields 14 reversals. 



Figs. 11-13 show the cycle counting results for the loading 
Cases 7-9, respectively. Cases 7-9 are asynchronous loading 
blocks, where the shear stress frequency is five times greater than 
the axial one. The difference between these three loading blocks is 
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based on the stress amplitude ratio (SAR) between shear and axial 
stress components. The block's SAR was 2, 0.5 and 0.6 in Cases 7-9, 
respectively. 



In loading Cases 8 and 9, the axial stress amplitude is great- 
er than the shear stress amplitude, i.e. axial damage predomi- 
nates over shear one. In Case 7 the opposite occurs, i.e. the 
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Fig. 10. Case 6: (a) von Mises stress space, (b) BS results with SWT, (c) BS results with F_Socie, (d) Block stress-time evolution, (e) WB results and (f) SSF time evolution. 



shear damage predominates over the axial one. In this case, the 
shear stress amplitude is twice the axial one. Cases 8 and 9 
have very similar results; therefore, little SAR variation causes 
little influence on reversals and accumulated damage results. 



Concerning Case 7 results, see Fig. 11, it can be concluded that 
SAR has a strong influence on the accumulated damage. How- 
ever, BS criterion yields 11 reversals in the three asynchronous 
loading cases. 
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Fig. 11. Case 7: (a) von Mises stress space, (b) BS results with SWT, (c) BS results with F_Socie, (d) Block stress-time evolution, (e) WB results and (f) SSF time evolution. 
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Fig. 12. Case 8: (a) von Mises stress space, (b) BS results with SWT, (c) BS results with F_Socie, (d) Block stress-time evolution, (e) WB results and (f) SSF time evolution. 



Concerning loading Cases 8 and 9, BS model extracts 11 
reversals in SWT and F-Socie critical plane criteria. Crack initia- 
tion estimations for Cases 8 and 9 were 24° and 26° in SWT and 



-20° and -19°, in F-Socie. Experimental results show -6° and 
-20° for crack initiation directions in Cases 8 and 9, respec- 
tively. Critical plane estimations without using a damage 
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accumulation criterion were, in Case 8, 0° in both SWT and F- 
Socie criteria. 



Moreover, in Case 9, the critical plane estimation without accu- 
mulative damage criterion were -23/+25 for SWT, and -20/+70 for 
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F-Socie criterion. The critical plane estimation in loading Case 7 
was 41° for SWT and 8° for the F-Socie criterion, respectively. 
The experimental crack initiation plane for loading Case 7 was 
-19°. Regarding Case 7, the critical plane direction without an 
accumulative damage approach were -39°; 38° in SWT and 
-13°; 12° in F-Socie. SSF reversals extraction from loading Cases 
7-9 were 2, 1.8 and 1.8, respectively. Rainflow method applied to 
the SSF time variation yields 19, 27 and 27 reversals for Cases 7- 
9, respectively. 

Loading Cases 10 and 11 results are summarized in Figs. 14 
and 15. These loading blocks are asynchronous loadings where 
the shear stress frequency is twice the axial one in Case 10 
and in Case 11 the axial stress frequency is twice the shear 
one. Concerning Case 10, BS criterion extracts 5 reversals in 
SWT and F-Socie at the most loaded plane (high accumulated 
damage). The critical plane estimation was -28° for SWT and 
15° for F-Socie. 

The Case 10 experimental crack direction was -18°. Critical 
plane estimations without an accumulative damage approach were 
±25 in SWT and ±17/±73 in F-Socie. WB yields 4 reversals and the 
most damaging reversal was identified in the last quarter of the 
loading period. Regarding the results of the Case 11, the BS yields 
the same reversals amount achieved in Case 10, i.e. 5 reversals at 
the most loaded plane. Critical plane estimations under the BS 
accumulative approach were -22° in SWT and ±17° in F-Socie cri- 
terion. In Case 1 1 the experimental crack initiation direction was 
-23°. The critical estimations without an accumulative approach 
were ±21° in SWT and ±23/±67 in F-Socie criterion. WB extracts 
4 reversals and the most damaging one occurs in the last quarter 
of the loading period. 

Regarding the SSF reversals extraction from loading Cases 10 
and 11 the results achieved were 1.8, and 1.5, respectively. Rain- 
flow results from the SSF time variation yields 14 and 23 reversals 
in the loading Cases 10 and 11, respectively. Table 2 summarizes 



Table 2 

BS critical plane estimations and experimental results. 



Case 


SWT (°) 


F-Socie (°) 


BS-SWT (°) 


BS-F-Socie (°) 


Experiment (°) 


1 


0 


±90; 0 


0 


0 


-5 


2 


±45 


±90; 0 


±45 


0 


0 


3 


±25 


±21; ±69 


±30 


±16 


29 


4 


±25 


±21; ±69 


±29 


±17 


31 


5 


0 


0 


4 


0 


0 


6 


-25; 23 


-23; 20 


24 


-20 


-22 


7 


-39; 38 


-13; 12 


41 


8 


-19 


8 


0 


0 


-2 


0 


-6 


9 


-23; 25 


-20; 70 


26 


-19 


-20 


10 


±25 


±17; ±73 


-28 


15 


-18 


11 


±21 


±23; ±67 


-22 


±17 


-23 



the critical plane estimations and critical plane experimental re- 
sults from loading Cases 1-11. 

Table 3 summarizes the results from the Bannantine and Socie, 
Wang and Brown, SSF and Rainflow approaches. The results are 
presented in reversals, two reversals are equal to one cycle. 

The BS cycle counting methodology aims to capture the damage 
behavior by using the rainflow method, a hysteresis loop method, 
to account block's reversals at each plane. 

In contrast, WB criterion yields more or less damaging reversals 
within a loading block without any connectivity to the physical 
damage process. 

It is not clear in the WB criterion why some regions of a loading 
block, i.e. reversals, are more damaging than other ones. However, 
BS and WB have similar results in the block's reversals extraction. 
The maximum difference between BS and WB reversals extraction 
was 3 in Case 8 and the minimum difference was 1 in several load- 
ing cases, see Table 3. SSF virtual cycle counting method in some 
loading cases extracts, from a loading block, a reversal quantity 
inferior to one, e.g. see Case 11 in Table 3. In Case 11 the SSF 
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Table 3 

Cycle counting results for Bannantine and Socie, Wang and Brown, SSF virtual cycle 
counting and rainflow methods, results given in reversals. 



Case 


BS-SWT 


BS-F-Socie 


WB 


SSF 


SSF-Rainflow 


1 


4 


3 


4 


3.4 


5 


2 


5 


3 


4 


3.7 


5 


3 


7 


7 


6 


4.1 


5 


4 


6 


7 


6 


4.0 


6 


5 


10 


7 


9 


6.9 


10 


6 


11 


11 


10 


1.9 


14 


7 


11 


11 


10 


2.0 


19 


8 


14 


11 


14 


6.8 


27 


9 


11 


11 


10 


1.8 


27 


10 


5 


5 


4 


1.8 


14 


11 


5 


3 


3 


1.4 


23 



method extracts 1.4 reversals, which is one reversal plus 40% from 
another one. To explain this result must be recalled the SSF virtual 
cycle counting paradigm, where the sum of the absolute values of 
peaks and valleys divided by two times the maximum peak yields 
the block's damage. The relative block's damage indicates how 
many times a loading block is more damaging than the reference 
one used to set up SN curves, e.g. sinusoidal R = — 1. Therefore, 
the SSF virtual cycle counting yields a non-integer reversal sum. 
The SSF cycles extraction amount was inferior to the ones obtained 
with BS and WB Brown criteria. 



15 

4.3. Fatigue life con-elation 

In order to estimate the block's fatigue life, it was used the 
reversals extracted in each cycle counting criterion coupled with 
a damage parameter and with the Miner's rule. 

Fig. 16 shows the fatigue life correlation between theoretical 
estimations and experimental results for BS criterion in association 
with a critical plane model. 

Sub-figures a and b show the BS-SWT results and sub-figures c 
and d show the BS- F-Socie results. BS-SWT criterion yields an 
accumulative damage parameter too low because the correlation 
tends to be above the boundary line. 

Regarding the F-Socie results, the correlation yields an incon- 
clusive trend in loading Cases 1-6 where can be found a big scatter. 
In contrast, for Cases 7-11, F-Socie fatigue life correlation yields 
very acceptable results with only two points outside of the bound- 
ary lines. 

The Wang and Brown fatigue life correlation is presented in 
Fig. 17. Sub-figures a and b show the fatigue life results for 
S=0.3, which sets up the Wang and Brown damage parameter 
to the infinite fatigue life estimations, and sub-figures c and d) 
shows the fatigue life correlation for half fatigue life estimations 
with S = 0.7. Likewise the BS results, loading Cases 1-6 show a 
big scatter for S = 0.3 and S = 0.7 and it can be seen that the accu- 
mulative damage parameter is too high in loading Cases 1-6. 
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Fig. 16. Fatigue life correlation for BS-SWT (a and b) and for BS-FSocie and (c and d) criteria. 
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Fig. 17. Fatigue life correlation for WB with S = 0.3 (a and b) and (c and d) with S = 0.7. 
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Fig. 18. (a and b) Fatigue life correlation for the SSF virtual cycle counting in association with the SSF equivalent shear stress. 



Nevertheless, for Cases 7-11, a better correlation with very Fig. 18 shows the proposed models fatigue correlation with 

acceptable results was achieved with S = 0.7, likewise observed experimental results where the scatter shown in BS and WB fatigue 
in BS criterion. life correlations at Cases 1-6 was reduced. 
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Moreover, the results for loading Cases 7-11 are within the fa- 
tigue life boundary with only 2 experimental results outside. 

Therefore, the proposed SSF virtual cycle counting method cou- 
pled with the new fatigue life criterion yields improved results in 
loading block's fatigue life estimations. 

Moreover, the new proposals have the advantage to be a very 
simple and quick method to estimate a loading block fatigue life. 
The complexity involved in the BS and WB reversals extraction 
leads to conclude that these cycle counting methodologies are 
not very practical to implement manually or even in a numeric 
code. 

5. Conclusions 

In this paper, one method and one criterion are proposed in or- 
der to characterize multiaxial fatigue under complex loading 
blocks. 

The proposed cycle counting method, the SSF virtual cycle 
counting, was implemented to extract the loading reversals within 
a multiaxial loading block and it is based in the SSF equivalent 
shear stress. 

The multiaxial criterion was implemented to estimate multiax- 
ial fatigue life under loading block conditions and it is based on the 
SSF virtual cycle counting method. In order to validate the pro- 
posed criterion it was carried out eleven loading blocks with differ- 
ent loading effects. 

Moreover, fatigue life correlations based in two multiaxial cycle 
counting methods, the Bannatine and Socie and the Wang and 
Brown, were compared with the new criterion fatigue life 
estimations. 

Results show that the proposed fatigue life criterion was suc- 
cessfully validated. Moreover, the proposed criterion has better fa- 
tigue life correlations than the ones achieved by the Bannatine and 
Socie and the Wang and Brown criteria. Furthermore, it should be 
underlined the simplicity to implement the new cycle counting 
method and the new fatigue life criterion. 
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